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Abstract. We present a description of CMB anisotropies generated by tensor perturbations in f{R) 
theories of gravity. The temperature power spectrum in the special case of f{R) = R" is computed 
^ using a modified version of CAMB package. 
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^ INTRODUCTION 

^ The cosmic microwave background (CMB) carries information from the last scattering 

^ surface that puts constraints on the multitude of proposed cosmological models and 

I— I the gravity theories they are based on. Amongst such theories are the f{R) theories 

^ of gravity which have become an interesting endeavour to correct for the degeneracies 

> of the concordance model [1, 2, 3, 4, 5]. 

(N 

^ THE PERTURBATION EQUATIONS 

Q In the 1 + 3 formalism [6], spacetime is sliced into constant time hyper-surfaces with re- 

spect to fundamental observers with 4-velocity u" = ^ where T is the proper time along 
the observers' world lines. Ua is time-like = — I). The metric tensor can be decom- 
posed into time and rest-space projection tensors: hab — gab + Ua^b- The derivatives in 
^ the 1 -I- 3 split are defined as: 

^ Projected covariant derivative : V aUJ- =1^^^ a^\- ■ -h' f ■ ■ dTe...^'" , (1) 

Time derivative : tb,,,''- = u^VaU...''- , (2) 

Generalized 3D curl : cm\Tab...c = u"ri^ne{a^'^Th,__cf, (3) 

where rjabcd is the totally antisymmetric tensor on spacetime. 

Kinematics are obtained by taking the irreducible decomposition of V^m^: 

VflMfc = ^Qhab + Oab + (Oab - Ua^b, (4) 
where = V'^Ua = 3H, Gab = O^^b) = ^ {aUb) = [h{a'' h)'^ - Ihabh""^] ^cUd, (Oab = 

= V[^M^], and Ab = Ub are the expansion scalar, the shear, the vorticity, and the 



acceleration respectively. H is the local Hubble parameter. A^, Oab, and (Oa characterize 
anisotropy in the expansion and therefore vanish in an exact FLRW universe [7]. 
The Riemann tensor, Rabcd^ can be decomposed into its trace and trace-free parts [8]: 



R 



(5) 



where Cahcd is the Weyl tensor and describes the vacuum contribution to the gravitational 
field. It can be decomposed further into the electric and magnetic Weyl tensors: 



^ab — W ^ ^acbd 



and Hab = I riacde u^u^Cb 



(6) 



The evolution of perturbations along the flow lines are described by seven propagation 
equations stated in [9], together with five spatial constraints. Here, we require only the 
linearised tensor mode equations characterized by: ^'^Eab — '^^Hab — '^^Oab — ^^^ab — 
0, leaving us with five propagations equations and one spatial constraint, of which we 
state: 

2 1 
= - 3 (^ab -Eab+^ ^ab and Hab = curl Oab . (7) 

Equations (7) determine E^b and Hab from the shear. Further differentiation and tensor 
harmonic decomposition of the shear equation in (7) gives: 



^-3(P + 3p) 



P%-^(P + 3p)0% 



(8) 



THE DYNAMICS OF f{R) COSMOLOGIES 

Metric variation of the f{R) and matter actions [10]: 

S= [ f{R)V^d''x + S„ 



'matter 



(9) 



leads to fourth order field equations (FE): 



1 1 1 

Rab - 7^ gab R = ^ Tab + -?r 
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^gabf+ (Va V^, - gabO)fR -\fR Sab R 



(10) 



where //? = ^ and Tab = -;^ ^gy is the energy momentum tensor (EMT) of 
standard matter. From the latter, the density, pressure, energy flux, and anisotropic stress 
are isolated via: 

_ Ttot,,a,,b „ ^^totiab „ Ttot,,ci„d _ rptot 

p = TabU U , p= -T^i, h , qa= -T^a U h a, T^ab = T^ab) (l^) 

qIii „m „m 

giving: p = ^+p«, p = ^+/, q^^'i^j^ql = ^ + ;r,^,. (12) 

JR JR JR JR 



The linearised form of p^, p^, q^, and 71^^^ are listed in [9]. We give here the most 
important contribution: 

^ab = Y [fj^R^{a^b)R - fRR OabR] (13) 
J R 

The f{R) Friedmann equations in flat FLRW with vanishing cosmological constant: 

lf-3fR- + 3-RfRR=Pm and \f - fR{H + 3H^) + -^(a^RfRR) ^ -Pm ^ 
2 a a 2 aat ^ ' 

CMB TENSOR ANISOTROPIES IN f{R) = R" GRAVITY 

For f{R) =R" and a single fluid with constant equation of state (EoS), the scale factor is 

given by: a(77) = 7^3(i+fi))-2n [10]. In the radiation dominated era: n"^ ^ ^0 and CO — ^. 
Separating the curvature and standard matter components, changing to conformal time, 
making the variable change = a'^Ok and choosing m = Eq. (8) becomes: 

4+[k^-2ri-^]uk^0 (15) 

for all values of n 7^ 2. This is the same equation of motion as for general relativity (GR). 

In addition, Equation (13) reduces to a direct relation between the shear and the 
curvature anisotropy: 

J, klfRR. k fn-l\R' 

ap /r \ P J R 

The evolution of the scale factor in i?" gravity may be obtained from Eqs. (14). 

RESULTS AND DISCUSSION 

Keeping a GR background allows to directly compare the influence of the first order 
modified equations. This is preliminarily our approach. Figure (1) shows the temperature 
power spectrum for f{R) = R", n = 1, 1.1, 1.2, 1.374, and 1.5. For comparison, the 
power spectrum for ACDM is also plotted. It was obtained by running CAMB [11] on 
the same background but with the original GR evolution equations. 

A few conclusions can be inferred directly from the power spectrum. We start by 
mentioning that the curve for GR and for i?" = /? are very close but not identical. 
The difference appeared after the use of the input file for the EoS. The plots for the 
two spectra are identical if a constant EoS for dark energy is used; typically CO = — I. 
Also, one notices that the features of the power spectra are shifted more and more 
towards small scales with increasing power of R. The departure from GR increases with 
increasing wavenumber and power of R as expected from the terms k and n — 1 in the 
numerator of the RHS of Eq. (16). The power difference between spectra for different 
values of n becomes larger after multipole ^ 160. Finally, the power decreases with 
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FIGURE 1. The temperature power spectrum for tensor perturbations in f{R) = R" gravity (n = 

1, 1.1, 1.2, 1.374, andl.5). The background is flat FLRW,Qfc = 0.035, ncDM = 0.315, Ha = 0.65, and 
-massive — 0- ^0 — 70km s ^Mpc ' . No secondary anisotropics or reionization effects were considered. 



increasing powers of R except for the interval 62 < / < 160 in the 77 power spectrum 

where the opposite happens. 

Calculation of the evolution of perturbations in R" gravity models with the exact 
cosmological scale factor is a work in progress. Equations (15) and (16) together with 
Figure (1) are the main result of this work. 
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